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psj , Abstract 

^ . By a quantum version of the Arnold transformation of classical mechanics, all quan- 

tum dynamical systems whose classical equations of motion are non-homogeneous linear 
. second-order ordinary differential equations, including systems with friction linear in ve- 

' locity, can be related to the quantum free-particle dynamical system. This transformation 

provides a basic (Heisenberg-Weyl) algebra of quantum operators, along with well-defined 
Hermitian operators which can be chosen as evolution-like observables and complete the 
entire Schrodinger algebra. It also proves to be very helpful in performing certain compu- 
tations quickly, to obtain, for example, wave functions and closed analytic expressions for 
time-evolution operators. 

a 

3 PACS: 03.65.-w, 02.20.-a, 2.30.Hq 
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(Sj '• 1 Introduction 

>. 

The description of the quantum damped harmonic oscillator by the Caldirola-Kanai model 
[TJ[2], which involves a time-dependent Hamiltonian, has attracted the attention of many authors, 
as could be considered one of the simplest and paradigmatic examples of dissipative system. In 
particular, its analysis from the symmetry point of view has proved to be very fruitful. In 
a purely classical context, the symmetries of the equation of the damped harmonic oscillator 
with time-dependent parameters were found in [3]. In [3] Cervero and Villarroel found, for the 
damped harmonic oscillator, finite-dimensional point symmetry groups for the corresponding 
Lagrangian (the un-extended Schrodinger group [5]) and the equations of motion (SX(3, M)) 
respectively, and an infinite contact one for the set of trajectories of the classical equation. They 
singled out a "non-conventional" Hamiltonian from those generators of the symmetry, recovering 
some results from [6J [7] . 

Some flaws has been associated with the quantum description of the damped harmonic 
oscillator according to Caldirola-Kanai equation. For instance, it is claimed that uncertainty 
relations are not preserved under time evolution and could eventually be violated [HJ|9], although 
this inconsistency seems to be associated with a confussion between canonical momentum and 
"physical" momentum |10| . 

There exists another interesting approach to the study of the classical damped harmonic 
oscillator, based on the observation that its classical equation of motion is a special case of 
the set of linear second-order ordinary differential equations (LSODE for short). In Classical 
Mechanics the family of solutions of a second-order differential equation corresponding to the 
motion of a given physical problem is sometimes related to that of a simpler system, considered 
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as a toy model, in order to import from it simple general properties which could be hidden in 
the real problem. Both physical systems should share global properties of the solution manifold, 
such as topology and symplectic structure. The paradigmatic example is the transformation 
described by Arnold in which brings any LSODE to the simplest form of the free Galilean 
particle equation. This transformation turns out to be extremely useful. In particular, it is 
possible to obtain the symmetry group of a particular instance of LSODE |3J, in which the 
symmetries of the action of the corresponding system can be found as a subgroup [12J. 

Therefore, it seems natural to try to generalize the Arnold transformation to the quantum 
level, to be denoted as Quantum Arnold Transformation (QAT), as much insight can be gained 
in the study of any system classically described by a LSODE and, in particular, the parametric 
oscillator or some of the systems which present dissipation. 

Several partial generalizations can be found in the literature. For example, in |13) . Takagi 
is able to provide a transformation which relates the Schrodinger equation of the harmonic 
oscillator to that of the free particle, and applies it to simplify the computation of the propagator 
by making use of the free one. |14) contains a slightly more general version (see formula (33) 
therein). |15] went a bit further considering the damped harmonic oscillator with constant 
parameters. None of them mentions the classical Arnold transformation, but it is underlying 
their reasoning. 

Implicitly, a generalization of the Arnold transformation was also contained in [7], the clas- 
sical version not being referred once more. It will be shown that some of their results formally 
converge with ours (see Section H~2j) . although they put emphasis on another aspects of the prob- 
lem, such as the analysis of unitarity and energy loss. Mostafazadeh [16J also pursued the idea 
of "connecting" different quantum physical systems by means of time-dependent unitary trans- 
formations, even representing arbitrary time-dependent diffeomorphisms |17] , His approach is 
rather general, but does not fully take advantage of the possibility of connecting with the free 
particle system and importing its symmetries. 

The relevance of generalizing the Arnold transformation is that it can be used to export 
properties from the free Schrodinger equation to that of the system corresponding to the given 
classical equation: the complete set of symmetries of the quantum free particle, the Schrodinger 
group [5j [18] , can be realized on the system under study, providing as many conserved quantities 
as in the free particle. In particular, the generators corresponding to the free position and 
momentum prove to be good basic quantum operators, constituting a quantization. In addition, 
the transformation turns out to be extremely useful to compute objects that would otherwise 
need laborious calculations, such as wave functions, the quantum propagator or the evolution 
operator. 

We shall see that time translations in the non-free system do not belong, in general, to the 
imported Schrodinger group. This is to be expected, as the energy in this system is not conserved 
when the classical equation of motion includes a friction term or a time-varying frequency. A 
deeper analysis of this fact can be found in [19] . 

To obtain the QAT, we begin in Section [2] by observing that a generalized version of the 
classical equation of the damped harmonic oscillator, in which the constant coefficients are 
promoted to be time-dependent, can be transformed by a local diffeomorphism of position and 
time, the Arnold transformation, into the equation of motion of the free particle. In Section [3] 
we generalize this interesting feature from the classical to the quantum theory. The Caldirola- 
Kanai equation for the damped harmonic oscillator is then a particular case of the general one. 
In Section U] we illustrate the use of QAT in a couple of simple examples presenting damping, 
the damped particle and the damped harmonic oscillator. 
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2 Classical Arnold transformation 




Mathematically speaking, the classical Arnold transformation converts any linear second- 
order ordinary differential equation (LSODE) into the free Galilean particle equation, that is, 
k = in 1+1 dimensions (we shall limit ourselves to this situation). 

From the physical point of view, the Arnold transformation relates the trajectories x(t), with 
initial conditions xq and p$ ~ xq, solutions of the LSODE, to those trajectories k(t) solutions of 
the free equation with initial conditions kq and ttq ~ Kq. Either (xo,J?o) or (^Oj^o) parametrize 
the common solution manifold M, and we shall adopt the unified notation (K,P). On this 
manifold, each physical system is characterized by the corresponding Hamiltonian as a function 
of K and P. The inverse of the corresponding Hamilton- Jacobi transformation then recovers 
the trajectories (x(t),x(t)) or (k(t),k(t)) out of the K,P variables. 

Following a similar notation to that in [3] , we give an overview of the Arnold transformation 
|llj . Firstly, let us recall that, given an arbitrary, non- homogeneous LSODE 

x + fx + uj 2 x = A, (1) 

where x = 4| and so on, and /, u) and A are arbitrary functions of time t, we can apply the 
transformation 

(2) 

Up being a particular solution of ([1]). We find that the differential equation above is transformed 
into 

x + f± + uj 2 x = 0, (3) 

i. e., every non-homogeneous problem is equivalent to a homogeneous one. 

The homogeneous Arnold transformation, is a local diffeomorphism which maps the free 
particle equation of motion into ([3]) : 

r T _ «!«) 

I U2 ^ , k = i — ► x + fx + u 2 x = 0, (4) 
[ K = 

where ui{t) and v,2(t) are independent solutions of Q. Applying the inverse diffeomorphism to 
the classical dynamical system Q, we can transform this equation into the free one. 

If we include external forces the transformation turns into the general Arnold transfor- 
mation, that we shall call simply A: 

( T _ m(t) 

_ , ^ = ^ x + fx + u 2 x = A. (5) 

l K ~ u 2 (t) 

This transformation could be understood as passing to coordinates analogous to co-moving 
spacial coordinate and proper time used in General Relativity, so that the system becomes "free", 
at least locally. 

Indeed this transformation is of local nature in time, in the sense that it is only valid for an 
open interval in time t. In fact, it can be shown that the equivalence in (j?]) and ([5|) is true up 

u 2 

to a factor -=§7, so that it holds in the interval where 112 does not vanish. This means that the 
transformation does not take the Euler-Lagrange operator associated with the LSODE itself to 
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that of the free system. For this reason it can not be claimed that both physical systems are 
actually equivalent. 

However, Arnold transformation can help to understand the physical system under study. 
In particular, as pointed in [3], it is possible to identify the set of contact symmetries for ([3]), 
and this way to arrive at the results found in [4J, which show the sets of symmetries for either 
the equations or the action from which such equations can be derived. 

It should be noted that, due to the general character of the transformation, we could miss 
the physical identity of position and time when performing such a transformation. But it will be 
possible to choose appropriate specific solutions u\ , U2 and u p with suitable initial conditions so 
that the identity of the variables is maintained^. We shall make use of this possibility in going 
to the quantum version of this transformation. 

Finally, having in mind the particular case of dissipative systems, we would like to remark 
that certain issues of the treatment of these systems are already apparent in the classical domain. 
For instance, time evolution is not a symplectomorphism, nor preserve the volume of the phase 
space. Obviously, the Hamiltonian function is not a Noether invariant. This becomes especially 
manifest and annoying when formulating the quantum theory. 

3 The Quantum Arnold transformation 

As already mentioned in the Introduction, several partial versions of the QAT can be found 
in the literature. Here we give a generalization that contains, as particular cases, those found in 

031 El EES Cl- 
in bringing Arnold's technique to the quantum world we must be aware, obviously, of the 
different philosophy of the quantum description and different nature of the equation of motion. 
The objects and structures that define a quantum system, namely the Hilbert space, the basic 
observables, the Hamiltonian operator, and the Schrodinger equation must be specified in a way 
that we are able to identify the same objects at both sides of the transformation. 

To this end, it is important to focus, in the free system, only on those operators corresponding 
to constants of motion, Noether invariants associated with its symmetry, that is, the Schrodinger 
group (which contains the centrally-extended Galilei group as a subgroup, containing in turn 
the Heisenberg-Weyl group of translations and non-relativistic boosts). This implies to fix the 
basic operators (that is to say, quantum operators which realize a unitary and irreducible repre- 
sentation of the common classical Poisson (Heisenberg-Weyl) algebra {K, P} = 1) so that they 
respect the Schrodinger equation, then having constant expectation values and being generators 
of the basic symmetrj|f| (the Heisenberg-Weyl sub-group of the centrally-extended Schrodinger 
group) . 

Those basic operators, in this form, are in principle the candidates to be related by a quantum 
version of the Arnold transformation, so that we shall have the situation as follows: On the one 
hand, a common Hilbert space Ti of wave functions ^(K) (L 2 (M)), which plays the role of 
initial values for both the solutions 4>(x, t) G %t of the Schrodinger equation relative to the 
quantum version of our original LSODE and those wave functions, ip(K,r) 6 Ti^, solutions of 
the free Schrodinger equation. And, on the other hand, the quantum Arnold transformation A 
relating Schrodinger equations and basic operators. As a crucial consequence, we shall obtain a 

For his restricted version of the quantum Arnold transformation, Takagi, in [13] . suggested that the transfor- 
mation could be specified by different solutions at different times, in principle avoiding the restriction of locality. 
However, care must be taken when making use of this freedom to prevent conflict with preserving the identity of 
the variables when it is desirable. 

2 By "basic symmetry" we understand, in general, those symmetries whose associated Noether invariants are 
enough to parametrize the classical solution manifold. 
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realization of the free symmetry on the quantum, non-free system. The following diagram can 
help to have a picture of the setup: 



U G (r) 



u(t) (6) 



1~Lq = 7~L y T~L = T~Lq 

i 

T-L^ (resp. Ht) is the Hilbert space of solutions of the free or Galilean (resp. non-free or 
corresponding to the LSODE) Schrodinger equation, Ug{t) (resp. U(t)) is the free (resp. non- 
free) evolution operator and 1 is the identity operator. Here, the Hilbert space % may be 
considered as the quantum analogue of the classical solution manifold, Ai, usually thought of as 
space of (classical) initial conditions. On ~H one must be able to measure all possibles physical 
observables in much the same way classical observables are characterized as real functions on 
A4, that is, functions whose arguments are constants along classical trajectories (functions of 
Noether invariants). 

The Hamiltonian of the non-free system, not being conserved in general, will not be related 
to any operator from the free particle. This is to be expected, since it is not a conserved quantity 
under the evolution of the physical system, neither at classical nor at the quantum level (in the 
sense that it does not have constant expectation values). It is important to remark that this 
implies that it is not possible to formulate a time-independent Schrodinger equation. 

More specifically, by extending properly the Arnold transformation (or the inverse) to the 
quantum case, we shall relate the space T~L^ of solutions of the free Schrodinger equation 

dip h 2 d 2 ip ,_. 

with corresponding classical equation 

K = 0, (8) 
to that space Tit where the quantum theory of the generic LSODE 

x + fx + uj 2 x = A, (9) 

is realized, the quantities /, cj and A being, in general, time-dependent. 

The classical equation ([9]) can be derived from a variational principle. We shall consider the 
Lagrangian 

L = e^-mi; 2 — —muj 2 x 2 + mAij (10) 

as our starting point. The Schrodinger equation can be derived from the corresponding classical 
Hamiltonian function 



H = — e~ f + (-mu?x 2 - mAx)e f , (11) 



according to the standard canonical prescriptions, leading to 



.d(f> h 2 _ f d 2 cj) ,1 2 2 a \ f , /x 

ih-— = e 1 — — =■ + ( — muj x —mAx)e J (p. (12) 

at 2m ax A 2 ' 
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For / linear in time, constant oj and A = this equation is commonly known as Caldirola-Kanai 
equation for the damped harmonic oscillator [TJ [2] . 

Even though both spaces of solutions of and (|12p . T~L^ and Tit respectively, will be related, 
and the basic quantum operators associated with the classical functions K, P realized as well- 
defined operators both on and TLt, we cannot still assure that both physical systems are 
actually equivalent. In fact, the evolution operator in time t, ifl-fjt, does definitely not leave 
invariant the space of the solutions of (fl~2|) in general, nor comes down to the space H, which 
means that it cannot be realized as an operator function of K, P (in sharp contrast to ih-§=, 
which is ~ P 2 ). We shall achieve the construction of well-defined Hermitian evolution- like 
generators imported from the free system via the inverse quantum Arnold transformation, but 
their eigenvalues, conserved indeed, do not correspond to the standard energjU. These operators 
close the Schrodinger algebra with the basic operators. 

Implicitly, this trick of considering operators different from the Hamiltonian to provide quan- 
tum numbers and obtain solutions of the Schrodinger equation as their eigenfunctions has been 
used extensively. For example, the operator found in [7] to label the energy-loss states, which 
coincides with the quantum operator H* corresponding to G5 in [H [20], turns out to be a gen- 
erator of the SL{2, R) subgroup of the Schrodinger group. We give here explicitly the frame in 
which this can be done: operators from the Schrodinger group can be chosen to play this role 
upon convenience. 

For the sake of simplicity, let us focus on the case with no external forces A = 0. The 
formulas corresponding to A 7^ are given in Appendix [Al 

The generalization of the classical Arnold transformation is obtained by completing (JJ]) with 
a change of the wave function. Explicitly, the quantum Arnold transformation, valid for every 
physical system with classical equation of the homogeneous LSODE type, is given by the (local) 
diff eomor phism : 

Ul(t) 

u 2 (t) 

^Kfi ( 13 ) 

i m 1 "2(0 -2 

(j)^u 2 {t)e 2 h w ^ »2W , 

where u\ and u 2 satisfy again the classical equation of motion in (x,t), ii\ = ^f,U2 = 
and W(t) = U\U2 — u\Ui = e.~^ ■ It is straightforward to check that by this transformation 
the Schrodinger equation of the free particle is transformed into (|12j) (with A = 0) up to a 
multiplicative factor which depends on the particular choice of the classical solutions u\ and 
v>2 (partial derivatives must be changed by the classical part of the transformation while wave 
functions are shifted by the quantum part). 

Now, we can impose on u\ and u 2 the condition that they preserve the identity of r and K, 
i.e., that (k,t) coincide with (x,t) at an initial point to, arbitrarily taken to be to = 0: 

«i(0) = 0, u 2 (0) = 1, «i(0) = 1, u 2 (0) = 0. (14) 

This fixes a unique form of the diffeomorphism for a given "target" physical system. However, 
the quantum Arnold transformation would still be valid if solutions u\ and 112 do not satisfy 
(|14p . The price to be paid would then be that the relation in the lower part of the diagram ([6]) 
above would no longer be the identity and basic position and momentum operators would then 
be mixed (see the end of this Section). 

3 The construction of a properly denned, Hermitian generator in standard time will take much effort requiring 
a sound analysis of the symmetry problem in damped systems. This has been done in |19| . 
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Formally, while the classical Arnold transformation in this case is: 

A: KxT — >RxT' 

( x ,i)_> (K,r)=A{(x,t)) (15) 

where T and T 1 are open intervals of the real line containing t = and r = 0, respectively, QAT 
can be written: 



A: Ht — ► 



(0 



0(s,t) i — ► v(«,r) = i(<£(x,t)) = ^(^(ije 2 » ww«aW* 0(x,t)) , 



(16) 



where .A* denotes the pullback operation corresponding to A. 

As already remarked, the basic symmetries of the free system are inherited by the LSODE- 
type system, as we are now able to transform the infinitesimal generators of translations (the 
Galilean momentum operator tt, corresponding to the classical conserved quantity 'momentum') 
and non-relativistic boosts (the position operator k, corresponding to the classical conserved 
quantity 'initial position'). They are, explicitly, 

* = ~> h i (17 » 

ih d . , 

K = K+—T—, 18 

m ok 

that is, those basic, canonically commuting operators with constant expectation values, that 
respect the solutions of the free Schrodinger equation, have constant matrix elements (and con- 
stant expectation values in particular) and fall down to well defined, time-independent operators 
in the Hilbert space of the free particle I? (R) . 

In general, these properties are satisfied whenever an operator 0(t) can be written as 

d(t) = u(t,to)du\t,t ), (19) 

where O is O(to) and U(t,to) is the evolution operator satisfying the Schrodinger equatiorjf]. 
If the Hamiltonian is time-independent, as in the free particle case, time-evolution is a one- 
parameter group and then U(t,tq) = U(t — tq). 

Defining a generic Schrodinger equation operator, S = ifr% — H, taking to = (U(t) = 
U(t,0)) and reminding that U(t,toyU(t,to) = 1, it is clear that, for operators of the form (|19|) : 

S 6{t) \4>{t)) = S U(t) 6\if>)=S U(t) \ip') = S = (20) 

l {x (t)\6(t)m)) = o (21) 

I (u(rf6(t)U(t)) = (22) 

±d(t) = ^ t d(t) + ^[H(t),d(t)} = (23) 

(where \ip(t)) = U(t)\i/j}), stating that those operators respect solutions, have constant matrix 
elements, fall down to define time-independent operators on the Hilbert space and are integers 
of the motion, respectively. (|1T|) and (|18p can be "de-evolved", so that 

U(rf ttU(t) = -ih^- , U(t^kU(t) = k, (24) 
Ok 



4 Note that 6{i) is not the usual Heisenberg picture version Oh (t) = U(t, to)^ 6 U(t, to) of its associated oper- 
ator in Schrodinger picture O, although their relation is very simple when the Hamiltonian is time-independent. 
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thus showing the properties above for the free particle. 

What we are doing is to focus on these integrals of motion, n and k, so that the new operators 
position X and momentum P acting on Ht are also integrals of motion in the non-free system. 
These will be the generators of the basic symmetry in the non-free system. Dodonov and Man'ko 
[7] obtained these operators in particular cases by direct calculation, imposing them to commute 
with the Schrodinger equation S. The difference is that, having related this system with that of 
the free particle, now it is clear how far one can go: the Schrodinger group and its enveloping 
algebra. Even more, the approach followed in [7] is only able to provide the basic operators 
(corresponding to linear functions in the classical solution manifold), since the condition of 
commuting with S, [S,0(t)] = 0, is more restrictive than that considered here of respecting 
solutions, which is equivalent to [S,0(t)] ~ S. 

Let us apply the QAT (|13p to (|17p and (|18p . For a given operator tt acting on Ti^ there 
is a corresponding operator P = A~ 1 -kA on Tit- The action on functions (t)(x,t) can then be 
obtained as follows: 

P(j){x,t) = A^TtA^t) = A- 1 nA*( y /uie 2 h w ^ <P{x,t)) = 

1 4_i2nJ_^3.^2 1 / Q , i m 1 u 2 2 

2 h w U2 A*~ l [-ih—A*(,/uie 2 ««2 (f>(x,t) 



U 2 v dn 



1 + 1 Hi J_ / Q im 1 "2 x 2 



2 ft W u 



U 2 

-ihu 2 ^--m^x) ( p(x,t) 



(fj i Tt% 1 Ul 2 2 \ 



(25) 



dx 

We can perform the same computation for the position operator and then we have: 

P = —ihu2T; mx-^- (26) 

ox W 

ill ih d , „s 

X = t^x + -ui— , 27 
W m ox 

thus providing the generators of the realization of the (centrally-extended) Heisenberg-Weyl 
symmetry on the physical system corresponding to a general LSODE. 

The properties of the operators (|26p and ([27p . i.e. preserving solutions of (fT2j) . having constant 
matrix elements, falling to the Hilbert space, are ensured by the properties of (|17|) and (|18p before 
the transformation, and will be explicitly checked for some particular cases in next Section. It 
will also become clear that the identity of both operators is preserved after the transformation. 

Apart from P and X , we can compute P 2 , X 2 and XP = \{XP + PX): 

P 2 = ~r?u\^ + ih 2 -^ x #- + m 2 §, x 2 + ih 7 ^ (28) 



■dx 2 W dx W 2 W 



W A mW ox m z Ox z mw 

+ n ft 2 d 2 uiu 2 + uiu 2 d u\ii 2 2 .,uiu 2 + uiu 2 , on . 

XP = —uiu 2 —-^ -ih — x- m— T x -ih — . (30) 

m ox 2 W ox W 2 2W 

Their first-order versions, valid on solutions of (11211 . are: 



2uiUi d 2 u\ d 2 .,u\Ui 



52 ..2mul d 2mu 2 u 2 d 2 ii 2 2 -u?u\ 2 mu 2 u 2 

p =lh -^di + th —^ x dx- + m —w 2 ~ x +lh ^r~ ( 31 ) 

X 2 = ^~j U " x 2 + l h 2 -^ x #- + ih^-i- + (32) 
W 2 mW ox mW ot mW 

vr> .,2mu 2 d uiu 2 + uiii 2 d iiiu 2 -uj 2 uiu 2 2 uiu 2 + uiii 2 

XP = —in — — in — x— m — ^ x —in — , 33 

W dt W dx W 2 2W ' v ; 
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which, together with X and P, close the whole Schrodinger Lie algebra: 



X,P 2 



P,P Z 



X 2 ,P 2 



2ihP 


mxp 



x,p 
x,x 2 
P,x 2 

X 2 ,XP 



ih 



-2ihX 
2ihX 2 



X,XP 
P,XP 
P 2 ,XP 



ihX 
-ihP 
-2ihP 2 . 



(34) 
(35) 
(36) 
(37) 



All these operators are well-defined on the solution space of the time-dependent Schrodinger 
equation, so that the action of one of them on a solution is again a solution. However, it is 
important to note once again that the Hamiltonian operator corresponding to the LSODE, that 
is, the quantum version of (|lip (A = 0), H, although being a second order differential operator, 
can not be expressed in terms of these operators in general and then it does not close a Lie 
algebra with them. But that which is worse, it is not even a well-defined operator on the space 
of solutions of the Schrodinger equation, T~Lf As a consequence, H is not the generator of a one 
parameter group corresponding to conventional time evolution. 

Instead, any linear combination of P 2 , X 2 and XP, say H , can be adopted as a well-defined, 
Hermitian evolution-like generator. It has an associated eigenvalue equation and real spectrum, 
and its eigenvalues can be used to label solutions of (|12p (A = 0) as its eigenfunctions. The 

particular choice of H to be taken depends purely on convenience and, for example, the similarity 
with the form of H of the particular physical system. 

We would like to point out that there is an essential difference between the approach followed 
in |20] and ours. The reason is that the un-extended Schrodinger group is considered there 
as the fundamental symmetry of the damped harmonic oscillator, the origin of which is the 
analysis of the classical equations of motion in |3j. The approach based on QAT provides 
directly a representation of a central extension of the Schrodinger group adapted to the specific 
LSODE-type systen@. For the relevance of central extensions in Quantum Mechanics, we refer 
to [Ml EH [221 [231 [21 [25]. 

Let us stress that QAT can be useful to quickly perform some calculations, avoiding tedious, 
direct evaluations which can become extremely involved in the system under study. For example, 
it can be used to compute the quantum propagator for any LSODE-type quantum system, 
following the idea of Takagi in |13j for the simple case of the harmonic oscillator, or even the 
evolution operator U(t), which becomes very difficult to evaluate exactly when the Hamiltonian 
is time-dependent and does not commute with itself at different times. 

Actually, the evolution operator of a LSODE system can be related with the free evolution 
operator. Having in mind the diagram Q, we write: 



(38) 



Here <p and (p are the same function of only one argument (k or x) and we will denote <p = (f> = ip. 
Then, 



U(t)iP(x) =A- 1 (U g (t)^(k)) 



1 i m 1 ^2 2 

e 2-w- 2 x A*~ l {U G (r)i}{K)) 



1 i m 1 &2 2 

_ c 2-w-^ A*- x {U g (t))A*- x {^(k)) . (39) 



5 This fact is of the greatest relevance for the analysis of the inclusion of time-symmetry in |19] for the damped 
harmonic oscillator. 
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To factorize the function tp and single out the general action of U(t), we compute 

.4*- 1 = ^) = e 1 ^ 1 ^)^! ^(x) , (40) 

where e 10 ^ 1 ^ 2 ^^ is a dilation operator which is not unitary. To unitarize this operator, the 
generator must be shifted from xj| to x-^ + ^, so that the true unitary operator is then 

tj n (-L) = etog(V«a)(*TC+5) = _J_ e los(V«a)*is (41) 
^v U2 y ^— v 7 

But the factor — 3= is already present in the previous expression of . Therefore, it now reads 



IU1 
i m 1 "2 „2 



= ___ e 27rw^ x e 2^ w i«2&?r e log(V«2)a;^ _ (-42) 

Its inverse is given by 

i » , x i /■ \ S ~ih„, „, d 2 im 1 «2„2 

Uit)- 1 = U(tf = ^e log(u2 ^e^ UlU2 e^e~^ w ^ . (43) 

Interestingly, we have been able to obtain an exact expression for the evolution operator as 
a product of operators. No perturbative approximation method, which could become cumber- 
some in some cases, is needed for any LSODE-related quantum system to obtain the evolution 
operator. These results hold for A^O (see Appendix [A|l . 

As a general comment before proceeding with the computation of the wave functions, let us 
go back to the relevance of conditions (fTi|) . Those have been chosen to preserve the identity 
of the variables and basic operators. Any other choice of solutions satisfying different initial 
conditions at any given initial time would have implications which must be kept under control. 
This was not taken into account in [7j, which might result in some confusing derivations. A 
general shift 

u\ —> au% + bv>2 , U2 — > cui + du% , (44) 

with the condition ad — be = 1 to preserve the value of the Wronskian W, is equivalent to the 
canonical transformation in the basic operators 

X -+ aX - —P , P -> -cmX + dP . (45) 
m 

That is, the freedom in the choice of the solutions u\ and it2, which is a SX(2,R) transformation, 
stands for a Sp(l,]R) transformation in basic operators and defines a family of quantum Arnold 
transformations. Then, the relation between Hilbert spaces in the lower part of the diagram ([6]) 
turns into a non-trivial transformation: 

Atfl){x) = lip{x) -»• v / de" ! ^#^(f ) . (46) 

Following the general ideas already noted, any linear combination of quadratic operators 
belonging to the realization of the Schrodinger group (any operator in the subalgebra of (|31D - 
(|33p ) can be chosen in such a way that its eigenfunctions solve (|12p . A specific combination of 
the operators (I3ip - ([33"1 ) with constant coefficients Cj and 7 was already chosen in [2oJf|: 

H* = —P 2 + -moj 2 X 2 + IXP . (47) 
2m 2 2 v ' 



6 For the damped harmonic oscillator with constant lo and 7, coinciding with uj and 7 resp., this operator is 
the only one from the 5*1/(2, R) Schrodinger subalgebra which commutes with the Hamiltonian. 
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The eigenfunctions of this operator, solutions of the Schrodinger equation, are 



faut) = i c 2h mx \ (« 2 - 7U1 /2)2+n2 U 2 + (ui-^mw i 



2mn . / 2m!i 



M 2 -7"i/2-tnm ) 2 (cd ( — V ft ) I C i? ( — ; J = 

V'(«2-1«i/2) 2 +!i 2 tt?' V 1 " W(«2-7«i/2) 2 +n2«2; 2 V v'(u 2 -7«l/2) 2 +A 2 «? 

(48) 

where C\ and C2 are arbitrary constants, D u are the parabolic cylinder functions |28j, T is the 



Gamma function, Cl = y oj 2 — and v is in general a complex number. 

In writing <f) v (x, t) we have kept the generality of the quantum Arnold transformation so that 
these solutions are valid for any LSODE-type system (the corresponding formula for a LSODE 
with a external force A 7^ is given in the Appendix [A]). This family of wave functions is more 
general than the one found by Dodonov and Man'ko in [7] in that it contains a set of functions 
valid when Z > Cj even for a general LSODE system. The associated spectrum of H* is 

h* = hCl(v + \). (49) 

To obtain these solutions, we have taken advantage of the QAT itself performing the fol- 
lowing steps. We solve the time-dependent Schrodinger equation corresponding to a harmonic 
oscillator with frequency Cl, considering both the attractive and the repulsive cases, so that we 
obtain solutions in terms of parabolic cylinder functions. Then, we take the QAT from this "in- 
termediate" system to the free one. We compose this QAT with the inverse QAT corresponding 
to passing the free system to the present LSODE system, obtaining this way solutions to the 
LSODE system Schrodinger equation, which are eigenfunctions of Hhou = 2m \ m Cl 2 X 2 . 
Finally, making use of the freedom in the choice in u\ and U2, we perform the shift 

Ml u 2 ->■ U2 — ^ui X— >X, P—t-P + m^X. (50) 
Its effect on the quadratic operators causes the expression of the particular combination 

H H on = + \m&X 2 -> 2^P 2 + \rrub 2 X 2 + |XP = H* (51) 

to change. As a consequence, the obtained solutions turn into (p u (x,t). 

The condition of normalizability must be imposed to retain the physically valid solutions. 
And we observe that the normalizability of the wave functions depends on the specific values of 
u and 7 in the expression (j47jll|. For u > 2, the normalizable solutions correspond to C2 = 
and v = n an integer. These functions are written then in terms of the Hermite polynomial:^. 
In the case when tl is imaginary, the solutions are Dirac-delta normalizable for v = — \ + iX, 
with A a real number. The operator H* shows a continuous, real, doubly degenerate spectrum 
in this case |29j (see also p 7 ] for constant 00 and 7 in the overdamping regime). The critical case 
f2 = can be obtained as a limit of this case. 

It must be emphasized then that the choice of these constants encodes the choice of the 
particular (arbitrary) quadratic operator belonging to the Schrodinger algebra used to label the 
solutions. In the framework of the quantum Arnold transformation, this freedom leads to other 
families of solutions, different from the one presented here. 



7 In the specific case of the damped harmonic oscillator, with constant uj and 7, it is possible to identify u) = lj 
and 7 = 7 (Cl = Q), so that the admissible wave functions vary depending on the regime. 

8 These are the eigenfunctions of the operator H* (there denoted as K(t)) found in [7] in the general LSODE 
case. It should be noted that they avoid, in this general case, the explicit mention of dimensional constants 
equivalent to Co and 7 and implicitly entrust the selection of the specific K(t) to the choice of the classical 
solutions, which might become rather confusing. 
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4 Dissipative systems: Hamiltonian vs. Hermitian operators 



Let us now have a close look at a couple of simple particular cases, extensively studied in the 
literature: the damped particle and the damped harmonic oscillator. Analogously, it is possible 
to analyze the harmonic oscillator from the QAT point of view. Although |13] and [7] contain 
some aspects of this analysis, it is possible to go a bit further and arrive at interesting results 

4.1 Damped particle 

For / = jt, a; = and A = in ()10p . a Lagrangian for the damped particle can be given: 

Ldp = ^me^x 2 , (52) 
where 7 is the damping constant. The equation of motion is then 

x + 7 x = 0. (53) 
Two independent solutions for this equation, satisfying initial conditions (|14|) , 

u 1 (t)= 1 ~ 6 , u 2 (t) = l, W(t) = e-*, (54) 
7 

provide the Arnold transformation for this system: 

_ l-e-T* 

7 

k = x (55) 
1<P = <j>, 

which turns out to be simply a reparametrization in time. 
The Schrodinger equation takes the form 

tK- = H DP <t> = -—e (56) 
and the corresponding basic symmetry generators 

P = -ih^-, jt = z + —(l-e*)£. (57) 
ox m.7 ox 

The crucial point is to realize that, in fact, the Hamiltonian operator H^p does not make 
sense as an operator acting on the space of solutions of (|56p . while P and X do. This can be 
checked by direct computation. For a given solution </>, the equation satisfied by </>' = Hop<f> is 
no longer (j56|) . the reason being that Hdp does not commute with the Schrodinger equation, 
while Xcj), for instance, does solve it: 

inl + ^-e-^—) (U) = (x+— (1-e-^)— + — e-^^) =0, (58) 
dt 2m dx 2 /V / V 7717 dx/\ dt 2m dx 2 J ' 

showing that the Schrodinger equation and X do commute. 

There is yet another way to check explicitly that Hdp is ill-defined in the quotient space 
by the time-evolution generated by itself. Formally, the equation (|56[) can be solved defining a 
time-evolution operator U(t,to). The fact that Hdp commutes at different times, 

[H DP (t 1 ),H D p(t 2 )]=0 (59) 
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makes the calculation of U(t,to) and its action on other operators simple: 
ttu 4- \ ~r fJ H DP (t')dt' « ( e -7*o- e -T*)^i 

U(t,to) = e h Jt o y ' = e 2m~,y >q£*, (60) 



If we choose to = 0, in agreement with conditions (|14|) imposed on solutions (|54jl . we recover the 
Arnold-transformed free evolution operator, as could be expected. However, the computation of 
the evolution operator directly and the possibility of obtaining it using QAT is not in any way 
trivial when f[59|) does not hold. 

By means of the action of U(t, 0) = U (t) on the basic operators P and X, or loosely speaking, 
using U(t) to "de-evolve" them until time to = 0, we can show that they match the form (|19jl 
and determine their action on wave functions depending only on x. This action, in this simple 
case, can be computed by expanding the exponential evolution operator and performing the 
commutation operations at each order of the expansion, leading to 

Q 

uUt)xU(t) = x, uUt)pU(t) = -ih— , (61) 

ox 

which do not depend on time t and take the usual Galilei forrrj^|. This automatically guarantees 
the three properties mentioned in the previous Section. In contrast, Hdp does not come down 
to the quotient by the time evolution generated by itself: 

UHt)H DP U(t) = ~^. (62) 

The reason for X and P to have good properties is that they are mapped from the free, 
basic symmetry generators (|18p and (|17p by the Arnold transformation (so that they are also 
symmetry generators of the damped particle system), while Hdp is n °t- It is then natural to 
map one of the quadratic operators belonging to the Schrodinger algebra of the free particle 
to make it act on the space of quantum solutions of the damped particle. This evolution-like 
operator defines a proper eigenvalue problem, that can be used to find solutions for (|56p . We 
could choose the operator H* already mentioned, but in this simple case we prefer to illustrate 
another rather natural possibility: the free Galilean Hamiltonian 

P 2 h 2 d 2 

G 2m 2m dx 2 ^ ^ 

In fact, the practical approach would be to solve the free, time-independent Schrodinger equation 
and Arnold-transform the solutions to obtain (non-stationary) solutions for (|56p . For example, 
plane waves are mapped into: 

(Pk{x,t) = e 2m ~i - ' . (64) 

The observations made above for the damped particle, being quite trivial, can help to clarify 
the general case. 



4.2 Damped harmonic oscillator 

Let us consider a friction function linear in time / = yt, as in the case of the damped 
particle, but also a non-zero constant frequency oj and no external force. The Lagrangian for 
the Caldirola-Kanai system reads 

L DHO = e^^mx 2 -^mu} 2 x 2 ) . (65) 

9 Note that the "de-evolved" operators take the form of those in (|57p when t = 0. But in general, the correct 
way to take the quotient by time evolution is that shown in ()6ip . 
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The classical equation of motion is then 



x + 7X + uj 2 x = . (66) 
Two independent solutions of (|66|) . satisfying convenient initial conditions (|14p are 

ux(t) = — e - 2*sinOt, u 2 (i) = e~ 2* cos Ot + -^e~2* s i n Ot, (67) 

where again W(t) = U\(t)u2{t) — U\{t)u2{t) = e -7 *, and 



o 

Note that these solutions have good limit in the case of critical damping u = \. 

Particularizing the quantum Arnold transformation (|13p for the free Schrodinger equation 
([?D, the Caldirola-Kanai equation is obtained: 

■U d< t> TJ A- ^ -ltd 2 <P 1 2 2 7 i - /« n 

th at =HDH0 ^ = -2^ e d^ + 2 muJxe +' m 

Basic quantum operators are now given by 

G 2 ^ 2 

— (20 cos fit + 7 sin fit) — + m— (7 2 + 40 2 ) sin fit x , (70) 

7* -yi 



20 v ' '&c 4fi 



e 2 e 2 d 

X = —— (20 cos Oi — 7 sin fit) x + i/i — — sin fit-— . (71) 
2fl y ' ' mfl dx y ' 

It is worth to note that these operators match those that were already found in [7] by hand, in 
looking for integrals of motion. 

Again, the key observation is that Hdho does not make sense as an operator acting on the 
space of solutions of (|69p . while P and X do respect solutions. Although this can be proved 
by direct calculation, it is more instructive to obtain the evolution operator U(t,to). But the 
fact that Hdho does not commute at different times [HDHo{ti)i Hdho^)] 7^ makes its 
calculation trickier using conventional methods, as already mentioned. 

Generically, one would approach the problem resorting to a perturbative method. An ap- 
propriate method to solve the operator equation for U (t) = U(t, 0) corresponding to (|69|) (we 
now make explicit the time dependence of the Hamiltonian) , 

ih^U(t) = H DHO (t)U(t), (72) 

is that of the Magnus expansion |27| (see Appendix IB"]) . 

However, it is a good idea to take advantage of the QAT instead. The explicit expression for 
the exact evolution operator encountered for the damped harmonic oscillator is rather involved 
and is found substituting ([67]) in (j4"2|) : 



U{t) 



2fie'2 I 



2fi cos fit +7 sin fit 



i m -2u) 2 e lt sin (It 2 
g 2 h 1(1 cos sin (It 

2 St 

^-i r e-T*sinfit(2ficosfit+7sinfit) 3- log(^, — 2 " e , 2 . n J x#- ,~n\ 

g 2m x ' dx^ g 2£2 cos Ut+7 sin fit } ox [ ( O) 
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W(t) is obtained analogously. It is then possible to check the de-evolution of the operators 
simply expanding the exponentials to the desired order, just as one would do for the free evolution 
operator. This allows to state that: 

(t) XU(t) = x, E/t (t) P U(t) = -ihS- , (74) 

ox 

as expected. 

The action of U(t) on HDHo{t) shows that it does not fall down to the quotient by the time 
evolution generated by itself. In fact, the de-evolution, 

uHt)H DHO U(t) = 

( h 2 d 2 1 2 2 \ / h 2 d 2 l 2 ,\ 

("2^ + 2 mUJ X ) ~ 7 H"2^^ - 2 mUJ X ) 

~/ 2 t 2 ( h 2 d 2 1 2 o 2u 2 , d , f .\ 3 . /r7rN 

+ V("^ + 2™ X + —^Tx ~ *)) + ° {t ) ' (75) 

does depend on time at any order in the time expansion. 

We remark once again that there is no actual need for a perturbative method to obtain U (t) 
in this case when QAT is used. 

The computation of the wave functions, solutions of (j69p . can follow the steps shown in 
Section [3l We select the operator H* particularized for the damped harmonic oscillator, that 
is, substituting u\ and 112 by (|67p and identifying Co = 00 and 7 = 7. Its general eigenfunctions 
are given by the corresponding expression (|48|l . and the spectrum (|49|) of H* will depend on the 
regime fixed by the specific value of il. 



5 Conclusions 

The analysis carried out in this paper permits to deal with the quantum theory of any 
LSODE-type dynamical system, using known properties of the quantum free particle. The 
quantum Arnold transformation provides basic operators and establishes that the symmetry 
group of the free particle, the Schrodinger group, can be transferred to a realization on the 
LSODE system. This result turns out to be of practical use when performing some computations, 
for instance finding solutions of the Schrodinger equation or the evolution operator, especially 
when the Hamiltonian does not commute with itself at different times. Even in these cases it 
is possible to give exact expressions, obtained in a non-perturbative manner. It is noteworthy 
that these calculations lead to the knowledge of objects in the quantum theory with the only 
requirement that the classical solutions of the LSODE are known. 

In a way, Arnold transformation allows to interpret LSODE-type forces, including dissipation 
linear in velocity, as effects observed in a "non-inertial reference frame". The use of the present 
scheme goes beyond the study of the simple damped harmonic oscillator, finding applicability 
in quite different branches of physics, such as Cosmology, where a scalar field appears (inflaton) 
satisfying equations in time which can be read as a LSODE. In this respect we are preparing 
the study of the specific example of the harmonic oscillator with time-dependent frequency |3Uj . 

It is worth mentioning that the Arnold transformation used in this work is just a particular 
case of a broader class of transformations which link free particle equations to even classical 
non- linear equations. Research in this direction would potentially lead to extremely useful and 
interesting results. We believe that a good starting point for this purpose was presented in |17| . 
although some effort to establish the explicit connection with the free particle would be in order. 
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Let us end up with a final general comment concerning the symmetry under time translation 
of the quantum system associated with a LSODE. Even though it is possible to set up a clear 
framework to deal with any LSODE-type quantum system by employing the quantum Arnold 
transformation, it does not provide by itself a well-defined operator associated with proper (true) 
time evolution. The reason is that conventional time evolution is not included in general in the 
symmetry group that can be imported from the free system: the Hamiltonian does not belong 
to the specific representation of the Schrodinger algebra. One may wonder what happens if 
time evolution symmetry is required. We pursue this interesting issue at least for the damped 
harmonic oscillator in |19] , 



A Appendix: Inhomogeneous LSODE 

We give here the general Arnold transformation with an extra external force term A and 
the corresponding generalization of the main results above. This computation follows analogous 
steps as those shown before. 

The general QAT is given by 



_ x-u p (t) 
K - u 2 {t) 



(76) 



The extra conditions to be imposed on the classical solution u p (t) to preserve the identity of x 
and t before and after the transformation are: 



u p (0) = 0, u p (0) = 0. 
In fact, the solution u p (t) can be expressed: 

u p (t) = K^mit) + K 2 {t)u 2 (t) , 
where: 



Ki(t) 



u 2 (t' 



^A(i')dt', 



/o W(t>) v ' ' ' w J W(t' 

The transformation (I76|) leads to the expressions for basic operators: 

= ~ lhU2 di ~ 171 w^ x ~ Up) ~ 171 w Up 

X = —Ax - u p ) + —u p H «i— . 

W W m ox 

The evolution operator reads 



(77) 



(78) 



(79) 



(80) 



U(t) 



%_ m 
,2 



J^ u , ur> a* _,. a_ 



(81) 
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Finally, the general solution of the corresponding Schrodinger equation, eigenfunction of the 
operator H* , is: 



4> v (x,t) 



^/v / 2^r(i/+l)^/(«2-7«i/2) 2 +n 2 uf 



(-Lrn(r-n \2 " 2 "1 /("2 -7"i /g) , t 2 "2~7"l/2 \ 



( irn f 1 („2.,(,\2 .;2\ M ""^"P ("2 / 2 ) , ("2 ~T"1 /2) amrrup \ 

e \2h J wv'VW " p ; ut fiVK(u2-7u 1 /2) 2ftW(u 2 -7"l/2) ^ W / 



U2— 7«i/2— irZui 
^/(«2-7«i/2) 2 +n 2 u 2 



H"7 



CiA, 



—f—(X-Up) 



^{u 2 -yu 1 /2) 2 +n 2 u 2 



+C 2j D_i_ 



. / 2m!) , > 

«Y — — (^-«p) 

^/(«2-7«i/2) 2 +n 2 « 2 



(82) 



B Appendix: The Magnus expansion 



The Magnus expansion was introduced as a tool to solve non-autonomous linear differential 
equations for linear operators and has the very attractive property of leading to approximate 
solutions which exhibit unitarity at any order of approximation. This is in contrast to the 
representation in terms of the time-ordering operator T introduced by Dyson. 

A solution to (|72p is given by 



U(t) = e {l( - t \ 0(0) = 0, 
and a series expansion for the matrix in the exponent 



(83) 



h(t) =J2& k {t) (84) 
fc=l 

which is called the Magnus expansion. We can write down the first three terms of that series: 
fii(f) = J* dt^-'-HDHoih)) (85) 

(86) 



'0 

1 rt 



fia(t) = o / dt i / dt 2 
1 Jo Jo 

]_ ft fh ft2 

fi 3 (i) = -J dti J dt 2 J dt 3 | 



jHDHo{tl), — jHDHofo) 



1 A r ^ a 2 * 

-rHr)Ho(ti) , --rHoHoih) , —tHdho{H) 
n In a 



+ 



^ r ^ a 2 

-jHr>Ho{h) , —jHDHofa) , —-^HDHo(ti) 



(87) 



However, a good iterative method to obtain the operator Cl(t) is given by the formula: 



Cl(t) = lim nWft), 

n— >oo 

where Ufc are the Bernoulli numbers and 

ad? (B) ee B , ad\(B) ee [A, B] , ad|(B) ee [ad^B), B] . 



(89) 



(90) 
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We have computed the operator 0(f) for the case of the damped harmonic oscillator to sixth 
order of approximation, to give: 

n^(t) = 

1* /7l 4. 7 2 * 2 , 7 4 * 4 ri , 2j£n , 7 6 * 6 fl , W , 32^ n W ft 2 S 2 1 2 2 \ 

- - f i + # + &(i + (-— A - -™ 2 A 

2 V ^ 12 t 360 v ^ 2m ^ 2 y 

+ + *d + m + «p + f# + w 1 )) - H (91) 

With this approximation to £l(t), obtained by means of (f88|) . one can obtain the quotient by 
time evolution of a certain operator 0(t), given by: 

fc=0 

These formulas lead to the same results found in Subsection 14.21 
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